ON SHAPE OPTIMIZATION AND THE POMPEIU 

PROBLEM 



ARUNAS GRIGELIONIS 

Abstract. The Pompeiu problem is considered as shape optimization 
problem. We show stability of the ball which is the minimum point 
of related domain functional. The proof is based on shape derivative 
method. Stability of the ball for general domain functionals invariant 
under the rigid motions is discussed. 



Introduction and statements 

Energy-type functionals F[u] = J J f (\x — y\) dco(x) du>(y), where u> is 
a measure of compact support, appears in statistical mechanics of systems, 
topological classification problems of knots, isoperimetric problems, har- 
monic analysis, discrete energy problems and other areas of pure and applied 
mathematics (see, for example, (HI E3 EE3 HI EE] )• 

Motivated by problems of shape optimal design HS1 EE3 El OH! we 
consider the model domain functional 

(0.1) F[n] = J J f(\x-y\) dxdy, 

n n 

defined on the set of bounded domains 0, C W 1 , n > 2, with smooth 
boundaries. Here / : (0, oo) — >■ M. is known function. 
Let V G C 1 (M n ,M n ) . Equation 

(0.2) x = V(x) 

will produce a flow : R x R n -> R n which moves a bounded smooth 
domain = S7o to its new position tijf = TY (J2 ) (see Fig. □ on p. inj- For 
all t close enough to zero the domain Q,Y nas been shown to be bounded 
and smooth |23j . 

Shape derivative in the direction of vector field V is defined by formula 
(0.3) Fp ;V]:=j t \ t=0 F[nY]. 
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Definition 0.1. We say nonempty domain SI is critical if 

F [SI; V} = 

for any V G C 1 (R n ,M n ). 

If SI is critical, then (see Section ^) equality 

(0.4) J f (\x — y\) dy = const 

n 

holds for all x G dSl with const = 0. 

As a special case equation (J0.4JI includes the inverse potential problem 

(0.5) J \x — y\ 2 ~ n dy = const. 

n 

This last equation defines the body which has constant gravity potential on 
its own shape. The observation based on the "moving planes" method was 
done in 

Proposition 0.2. Let C 1 function f is strictly monotone. If nonempty C 1 
domain SI solves equation \0.J\) , then it is a ball. 

Thus, if C 1 function / is strictly monotone, then any critical point of 
the domain functional is the ball of fixed radius. If function / is not 
strictly monotone, then critical domains of different shapes are possible. 

Usually stable minimum point of general domain functional F is defined 
as the domain SIq at which inequality F [Slo] < F [SI] holds for all SI close 
enough to SIq (see, for example, jHj). Even despite of the difficulty in proper 
understanding of the closeness, domain functional 1)0. lj) has no stable mini- 
mum (and maximum) points in this sense because it is invariant under the 
rigid motions. 

To reach the uniformity with definition of critical point we shall introduce 
the notion of weakly stable minimum point for the domain functional F. 

Definition 0.3. We say the domain SI = SIq is stable minimum point of the 
domain functional (|0.1|) in the direction of vector field V £ C 1 (IR n ,IR n ), if 
either for all t, t ^ and close enough to zero, inequality F [Sl ] < F [fif] 
holds, 

either for all t close enough to zero, SIq as the rigid body coincides with 



The domain Si = SIq is said to be weakly stable minimum point of the 
domain functional IJ0.1|) in C 1 (R ra ,lR n ), if this domain is stable minimum 
point of F in the direction of vector field V G C 1 (R n ,lR n ) for any vector 
field V G C 1 (M n ,M n ). 
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Remark 0.4. a) Inequality 

(0.6) F[Q ; V, V] > 0, 

here F denotes the second shape derivative of the domain functional F, is 
sufficient for the critical point 0, = f?o to be stable minimum in the direction 
of vector field V. 

b) If V = const, then any Qo as the rigid body coincides with Q}[ for all 
t £ M. 

c) If normal part of V on the boundary dflo is zero, then = f° r an 
t close to zero. 

Statements b) and c) are easy to understand from geometrical point 
of view. They follow from the standard properties of ordinary differential 
equations. 

If C 1 function / is strictly monotone, then by Proposition 10.21 the domain 
functional F has at most one critical point to be considered as a rigid body 
which, of course, is weakly stable in C 1 (M n ,M n ). 

In this paper we consider domain functional (|0.1[) when function / ( | • | ) is 
positive definite, that is 

(0.7) ff f{\x-y\)u(x)-u(y) dxdy>0 

for all complex u € C§° (W 1 ) and thus inequality 
(0.8) F [Q] > 

holds for all f2 bounded. 

Theorem 0.5. Let F is the domain functional \U. 1\) generated by positive 
definite function f (|-|) and 

(0.9) F [O] = 

for some f2 the ball. 

Then 17 is weakly stable minimum point of F in C 1 (R n ,M n ). 

Investigation of the domain functional (|0,1[) with the positive definite 
function / (|-|) is stimulated because of its relation to the Pompeiu problem 
[221 • We wil1 P rove 

Proposition 0.6. Let £1 be bounded nonempty domain. 

a) Iftt satisfies hU. y\) for some positive definite function / (|-|), then there 
is A > 0: 

(0.10) in (0 = for all £ € M. n , |£| = A, 

here xn denotes the Fourier transform of indicator function xn- 

b) On the other hand, if condition I0.1fl\) holds, then F [CI] = for some 
positive definite function f (|-|). 
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Domain Q C M n which satisfies ()0.10f) is known as domain without the 
Pompeiu property 5 . If its boundary is smooth and connected, then 
overdetermined problem 



has nontrivial solution [21110]. It is expected (Schiffer conjecture [2T]) such 
the domain to be necessarily a ball. 

In ^ ^] the Schiffer conjecture is affirmed for "small" perturbations 
of the ball. Kobayashi ^3] use precise Fourier transform estimates of do- 
main indicator function. Agranovsky and Semenov PP, instead, exploit the 
overdetermined problem (|U.11|) . Their approach leads to the generalization 
of results in Riemannian spaces. More on the Pompeiu and Schiffer problems 
see the expository paper 

The proof of the Theorem 10. 51 mav be derived from [IJIE1- Instead of, we 
consider the Pompeiu problem as shape optimization problem and shall give 
the direct proof of the Theorem 10.51 based on the shape derivative method. 

The paper is organized as follows. In Section ^ shape derivatives of the 
domain functional (|(J.1|) are calculated. The domain equation (|0.4[) which 
describes critical domains is derived and an expression of the second shape 
derivative of the domain functional (|0.1j) in a case of positive definite function 
/ (|-|) is examined. 

The proof of the Theorem 10.51 is preceded by 

Lemma 0.7. Let F is the domain functional \0.1]) generated by positive 
definite function /(|-|) and F = for O, the ball centered at x = 0. 
Then F [Q; V, V] = for vector field V if and only if the corresponding 
function v = (V,x) on the boundary of £1 coincides with first order spherical 
harmonic. 

This lemma is proved in Section [2j 

Based on the Lemma 10.71 the proof of Theorem 10.51 is done in Section |31 
and consists of two steps. First, we consider special choice of vector fields 
V when all domains QY has the same center of mass x = (Lemma 13.1(1 . 
Second, the case when center of mass moves, is reduced to the previous step 
by suitable choice of non-autonomous vector field. 

The proof of the Proposition 10.61 about the Pompeiu problem equivalence 
to the shape optimization problem is done in Section 0J 

In Section |S] the remark on generalization of Theorem l0.5l to other domain 
functionals is done. Then a few numerical examples of shape evolution to 
the solution of the Pompeiu problem are presented. 



(0.11a) 
(0.11b) 
(0.11c) 



Au + X 2 



u = in O 



u = 1 on dfl 
= on dQ 
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Figure 1. 

1. Shape derivatives of the domain functional 

Let / : (0, oo) — ► R be the smooth function. We consider the domain 
functional 



(1.1) 



F[n} = J j f{\x-y\) dxdy. 



n n 



Then, according to notions on page ^ (see also Fig. [J), 
(1.2) 



F[nY]= J J f{\T t v (x)- T t v (y)\) det DT t v (x) det DT t v (y) dxdy, 



here det DT? denotes the Jacobian of the mapping T% 
considered as parameter. 
Using equalities 



pn . won 



(1.3) 
(1.4) 
(1.5) 



dt 



±T? (x) = V (T? (x)) , 
T? ( X ) = x , 
\ t=0 det DT? (x) = div V (x) , 



(see [231 on calculus technique) and verifying the equality 



(1.6) T At= f(\T?(x)-T?(y)\) 



dt 

one gets 
(1.7) 



(V,/ (|x — y\),V (x)) + (V,/ (\x — y\),V (y)) 



-\ t=0 F[nr]= j IW*fQx-v\),V(x)) + (V v fQx-y\),V(v))] 



-f (\x - y\) [div V (x) + div V (y)]}dx dy 



div x [f(\x-y\)V(x)} dxdy = 2 



dSlo Lfio 



f(\x-y\) dy 



(y,7]) (x) ds x , 
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here rj is exterior unit normal vector field on d£lo. 
If gives an extremum to F, then necessarily 



(li 



/(|x-j/|) dy = 



n 



for all x E dtlo. 

Replacement of the "initial moment" t = by the arbitrary moment t 
leads to equality 



an} 



f (\ x - y\) dy 



(V, V [)(x) ds x , 



here j]( denotes exterior unit normal vector field on dQY ■ 
Now let us calculate the second derivative 

(1-10) F[n ;V,V]:=-^\ t=0 F[nY] 

in case when f^o is critical domain to functional F. 

Applying formula (5.9) in page 1105 of and because of 1)1.8(1 one gets 

(1.11) 



ion} 



f (\ x - v\) dy 



(V,TI?)(X) ds x 



J {j t \ t=0 J f{\T^(x)-T^(y)\)detDT^(y)dy 



(v,vY) (tY{x)) 



+ (n-l)K(x) J f(\x-y\) dy ■ (V,r)) (x)}ds x = 
[div y {f (\x - y\) V (y)) + {V x f (\x - y\) , V (x)>] dy \ (V, r,) (x) ds x = 



f(\x-y\)(V,ri) (x) (V lV ) (y)ds x ds y 



+ J (V x Jf(\x-y\)dy,V(x)UV,r))(x)ds x , 

here K [x) denotes mean curvature of 8Qq at the point x. 

If function / (|-|) is positive definite and F [Qo] = 0, then indicator func- 
tion XQ niinimizes functional 

(1.12) F[<P]=J Jf(\x-y\)<p(x)<p(y)dxdy 
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defined on the set of compactly supported functions. Then necessarily 

(1.13) J f(\x-y\)dy = 

n 



for all x G R™. As a consequence, the second term in the last part of (|1.11|! 
is zero. 

So i£F[Q ] = 0, then 



(1.14) F[n ;V,V} = 2 J J f(\x-y\)v(x)v(y)ds x ds y , 

for all V G C 1 (M n ,M n ) . Here 

(1.15) v (x) = (V, rf) (x) 
denotes the normal component of vector field V on <9f2o- 

2. Proof of Lemma IUT7I 

The proof is based on standard properties of Fourier transform and of 
Bessel functions. For them we refer to ^5] and |19| . 

Let B be a ball of radius R centered at x = and F [B] = for the 
domain functional (jU.lj) with positive definite function / ( | • | ) . Then 

(2.1) 

oo 

= 1 1 f(\x-y\) dxdy = I J \xb\ 2 (0 d^d/x(r) = 
B B |£|=r 

(2nR) n J J 2 ^ |ra ds ? d/i (r) = (27ri?)" oj n ^ J r' l jl (Rr) dfi (r) , 
|£|=r 

here io n -\ denotes the surface measure of the unit sphere. 

Jn (t) 

Because function \ is entire and fi is positive Borel measure of poly- 
* T __ 

nomial growth |16j . equality in ()2.1|) is possible only if measure fi is discrete 
and 

(2.2) supp [i C |t|J« (Rt) = o} . 

Now, let us consider equality F [B; V, V] = 0. We have 

F[B;V,V}= [ [f (\x - y\) v (x) v (y) ds x ds. 

(2.3) 



y 

8B 8B 



\(XdBv) A \ 2 (0 ds c dv(r). 



o l€l= 
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Function v £ L 2 (dB) enables decomposition 

00 

(2-4) v(R0 = ^2c k Y k (0, |f| = l, 

fc=0 

for spherical harmonic Y k of order k [1 
Then 

00 



F[B;V,V}= J r n -< J ]T \c k \ 2 \ ( XdB Y k ) A \ 2 (£r) da (r) = 
= j rn ~ l j El^l 2 ^^ Y k \Rr) ds^{r) 



|f|=l k=0 

since 

/J n + 2fc-2 (-Rr) 
n (0 e^'*) d^ = const ^ y fc (x) 

|C|=i 

and cons/; 7^ jTH]. 

Bessel functions of different order has nonintersecting sets of zeroes |19j . 
so due to relation ()2.2j) equality 



(2.7) J J f(\x-y\)v(x)v(y) ds x ds y = 

is valid if and only if 

00 

f (\x - y\) v (x) v (y) ds x ds y = |ai| 2 w n _i / rjl(Rr) d/x (r) 



(2.8) 9B9B 

y y f(\x-y\)Y l (x)Y 1 (y) ds x ds y , 
dB dB 

where Y\ is the first order spherical harmonic. 
Remark 2.1. Function Y\ has the form 

(2.9) Yy (x) = b\X\ + b 2 x 2 + ... + = (b, x) , 

i.e. on the surface of the ball B it coincides with a normal component of 
a constant vector field. 
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3. Stability of the ball 
In this Section we shall prove Theorem 10.51 

Because of Lemma 10 .71 Remarks 10 . 41 and 12 . 1 1 it is sufficient to verify stabil- 
ity of the ball f2 = centered at x = in direction of vector field V when 
function v (x) = (V (x) , x) on the boundary dQ is the first order spherical 
harmonic. 

Lemma 3.1. Suppose that point x = is mass center of the domain Q.Y f or 
all t close enough to zero and function v (x) = (V (x) , x) on the boundary 
of the ball Qq is the first order spherical harmonic . Then v = on dtlo . 

Proof. For all t close enough to zero and for any i = l,n equalities 

Xj dx = 



and 



d 

|i=o 



x i^ x = J Xiv(x) ds 

holds. Consequently, v = on the boundary 80,q. □ 

In general case, mass center xt of the domain flY is not zero. We consider 
the flow T t (x) := T t (x) — xt- This flow generates the family of domains 
Q.Y which coincide as rigid bodies with fij 7 and are centered at x = 0. Of 
course, the new vector field 

V(t,x) :=^Zi V (s) 

is non-autonomous. Fortunately, shape derivative method is applied also in 
the case of non-autonomous vector fields V G C° (/, C 1 (R n ,R n )), I C R, 
G I |23j . Moreover, withdrawal of formula (jl,14|) in Section ^remains the 

same in the case of v (x) = (V(0,x),rj (x) 



4. Shape optimization and the Pompeiu problem 
In this Section we shall prove Proposition 10.61 

Proof of a) Let equation (|0.9|) has bounded domain O as its own so- 
lution for some / ^ 0. By Bochner-Schwartz theorem |16j and because of 
spherical symmetry of f(|-|) , there exists [i — positive Borel measure on 
(0, oo) of polynomial growth: 

F[Q} = J Jf(\x-y\)dxdy = J \xn\ 2 (0 d/i(|£|) = 

oo 

J J \xn\ 2 (0 ds 5 d^(r) , 



q n 



(4.1) 

|e|=r 
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Xn(0= / e^dx. 



where 
(4.2) 



Then equality 

(4.3) F [O] = 

implies the existence of A > : 



(4.4) 



J \M 2 (Z) ds § = 0. 



IC]=A 



Consequently, 

(4.5) xn (0 = 

for all |£| = A. Because of xn (0) = mes ^ > 0, one has A > 0. 

Proof of b) Suppose that O has no the Pompeiu property. Then (j2j, 
0) there exists A > equality (|4.5() is true and, therefore, 

(4.6) 







\xn\ (0 dst 



a n 



l€l=A 



l€l=A 

dxdy = (2vrA) : 



n n 



e i{i,x-v) dxdy 
J n-2 (A \x — y\) 

2 

n — 2 

\x — y\ 2 



dsc 



dx dy, 



here J p denotes the Bessel function. 

This will end the proof of Proposition IU.61 



5. Concluding remarks 

1. Let F be a domain functional defined on the set of bounded domains 
O C W 1 and invariant under the rigid motions of W 1 . Suppose that for the 
ball B centered at x = the condition 

(*) if F [B; V, V] = 0, then function v (x) = {V (x) ,x) is first order 
spherical harmonic 

holds. 

Proposition 5.1. IfF[B; V, V] > for all V G C 1 (M n ,IR n ) and the con- 
dition (*) holds, then the ball B is weakly stable minimum point of domain 
functional F in sense of Definition \U.°<A 



2. Concerning the Schiffer problem. The result of Theorem 10.51 does 
not cover ones stated in I13j . This is because spectral parameter A in 
(|0.11[) may vary with the domain 0. Though this more general situation 
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Figure 2. 

may be considered on the basis of Lemma l().7l the proof requires additional 
techniques. 

3. Methods based on the shape derivative allow numerical experiments 
in the Pompeiu problem. In Fig. |2] evolution of the long thin ellipse and 
of the quadrate to the "figure" without the Pompeiu property in the shape 
antigradient direction is presented. 
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